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Abstract 

We consider a massive relativistic particle in the background of a gravitational 
plane wave. The corresponding Green functions for both spinless and spin 1/2 cases, 
£>. , previously computed by A. Barducci and R. Giachetti [1], are reobtained here by 

IQ j alternative methods, as for example, the Fock-Schwinger proper-time method and 

the algebraic method. In analogy to the electromagnetic case, we show that for a 
gravitational plane wave background a semiclassical approach is also sufficient to 
provide the exact result, though the lagrangian involved is far from being a quadratic 

r — ' 

q ■ one. 

S3 
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$*• 1 Introduction 

si ' 

Green functions are basic ingredients in quantum theories. Particularly, they are of great 
importance in the computation of scattering amplitudes, as well as atomic energy levels. 
However, only a few exact results for Green functions of relativistic particles or even non- 
relativistic ones under the influence of external fields are available in the literature. Of 
particular importance, among the methods of obtaining Green functions is the so called 
Fock-Schwinger method. It was introducted in the context of relativistic quantum field 
in 1951 by Schwinger [2]. It has since been employed mainly in relativistic problems 
such as the calculation of bosonic [3] and fermionic [4, 5, 6, 7] Green functions in external 
fields. The most common exact solutions are those involving the non-relativistic Coulomb 
potential [9], the relativistic Dirac-Coulomb potential [10], a constant and uniform elec- 
tromagnetic field as well as the electromagnetic field of a linearly polarized plane wave 
[2, 11] or even particular combinations of constant fields and plane wave fields [12, 3]. 

In the presence of an external gravitational background the problem becomes ex- 
tremely complicated and only a few solutions using algebraic techniques or mode summa- 
tion method are known. In this context, a simplifying assumption is to consider a weak 
gravitational field. Using path integral and external source methods it is possible to cal- 
culate the Feynman propagator in the case when the gravitational field is that of a plane 
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wave[13, 14]. Recently, Barducci and Giachetti [1] have considered the wave equations 
for spin and spin \ particles in a weak external plane wave gravitational field. They 
obtained for the wave functions Volkov type solutions [16]. A similar ansatz for the Green 
functions is verified to be correct. These results have a resemblance to those for a charged 
particle in a plane wave external electromagnetic field. Since the latter case has been 
treated in a deductive manner by the proper-time method [2], it is interesting to apply 
this technique for the case of a weak background of a gravitational plane wave. 

The purpose of this paper is to provide many alternative methods for computing the 
above mentioned Green functions. This work is organized as follows: in section 2 we obtain 
the Green function for a spin massive particle in the presence of a weak background of a 
gravitational plane wave by using the Fock-Schwinger method. In section 3 we apply this 
same method and obtain the corresponding Green function for a spin 1/2 particle. Then, 
in section 4, we show that these solutions may be constructed in a much simpler manner 
by an operator technique. In section 5 we reobtain the bosonic Green function by trying 
a convenient ansatz. In section 6, we show that a path integral semiclassical aproach is 
sufficient to yield the exact results, even though the lagrangian involved is far from being 
a quadratic one. Finally, section 7 is left for conclusions and final remarks. 



2 Proper-time method for a scalar particle 

In the linear approximation the gravitational field is described by a small perturbation to 
the flat metric, namely, 

g ^(x) = rT-h^ix) (1) 

where r/^ is the flat space metric and is a small perturbation. The imposition of 
harmonic condition in the linear approximation gives 

d^ v {x) = l -d v h{x) . (2) 

It is convenient to choose a restriction of the harmonic gauge and impose the conditions 

d^ v = , h(x) = . (3) 

Finally, the linearized gravitational field is taken to be one produced by a wave of arbitrary 
spectral composition and polarization properties, but propagating in a fixed direction so 
that 

Ku{x) = a^F(0 (4) 

where £ = n-x and the propagation vector is a light-like one, which satisfies the condition 
n 2 = 0, and F is an arbitrary function. The conditions imposed on give 

nX v = , Tra = . (5) 

We note that linearization implies that terms involving second and higher powers of a are 
to be dropped. In the following calculations it may be convenient sometimes not to do 
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this in the intermediate stages. With these assumptions, the Green function for a scalar 
particle in a weak external linearized gravitational field satisfies the equation [13] 

_ hrdpdv + m 2 ]G(x -y) = -5(x - y), (6) 

where h^ v (x) is of the form described previously. In order to apply the proper-time method 
we write the Green function in the form 

roo „ 

G(x', x") = -i ds e~ im s (x'\e- isH \x") , (7) 
Jo 

where the proper-time hamiltonian is given by 

H = -p 2 +p»a, u p»F(0, (8) 

with = id^ and 

\tf,tf\ = -iif. (9) 

We adopt Schiwnger's notation, in which unprimed quantities are used for operators 
while primed quantities are used for the corresponding eigenvalues. In this sense, the 
eigenvalue equation for the operator x^ is written as 

x »\ x ') = x >»\ x i) , (10) 

where, for simplicity, we omitted indices in the eigenvector \x'). Next, given an operator 
O, we introduce the operator O(s) defined as 

0(s) = e lsH Oe- isH , (11) 

which satisfies the Heisenberg-like equation of motion 

i±0(s) = [0(s),H\. (12) 

Of course, O(0) = O, so that p^, etc., mean the same as p M (0), and x M (0) and so on. 
Hence, we have 

x^s) (e isH \x')) = x' (e isH \x')) , (13) 
so that we conveniently define 

\ x 's) = e isIi \x') . (14) 
Note that from (15) and (11) we also have 

[x< i (s),jr(s)] = -irf v - (15) 

Using this "Heinsenberg picture" , the Green function can be cast into the form 

roo . 

G(x',x") =-i ds e~ im s (x's\x"0) , (16) 
Jo 

where the "Schrodinger-like" propagator (x's\x"0) satisfies the differential equation 

id s (x's\x"0) = (x's\H\x"0) , (17) 
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submitted to the initial condition 



\im(x's\x"0) = 5(x' -x"). (18) 

Schwinger's method for computing (x's\x"0) basically consists in the following steps: 
(%) we first solve the Heisenberg equations for operators p fl (s) and rr M (s) and write the 
proper-time hamiltonian in terms of x^(s) and 2^(0), instead of 2^(0) p M (0); (%%) then, us- 
ing the commutator between x^(s) and 3^(0) we write conveniently this hamiltonian 
in a ordered form in time s, namely, with operators x^(s) in all terms put on the 
left side, so that equation (17) can be imediately integrated in s to yield (x's\x"0) = 
C(x',x")exp{-iJ s F(x',x";s')ds'}, where F(x',x";s') := (x's\H\x"0) / (x's\x"0); (Hi) fi- 
nally, the integration constant C(x',x") is obtained by imposing the constraints 

(x's\p^s)\x"0) = i^-(x's\x"0) , 

OX ^ 

(x's\p tl (0)\x"0) = -i-— (x's\x"0), (19) 

ox" 

as well as the initial condition (18). A pedagogical introduction of Schwinger's method 
can be found in [17, 18] (for other applications to non-relativistic problems see references 
[19, 20]) 1 . The Heisenberg equations of motion are given by 

j-x^s) = 2[r 1 \-a» v F(i{s))]p»{s) (20) 
j- s jf{s) = n^F'(as))p a (s)a a ^(s), (21) 

where in the last equation F' means dF/d^. In order to solve previous equations for x^(s) 
and observe initially that they imply the following ones 

±Z(s) = 2n^(s)=:2n-p(s), (22) 

^(n-p(sj) = 0, (23) 

^-(p»{s)a^p v {s)) = 0. (24) 
Equations (22) and (23) lead to 

»-PW = ^a (25) 

and equation (24) allows us to write 

p^s)a^(s) = C, (26) 



1 Schwinger's method for the non-relativistic oscillator was developed independently by M. Goldberger 
and M. Gellmann in 1951 in the context of statistical mechanics [21]. 
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where C is an operator that does not depend on s. In order to integrate equation (21), 
we multiply and divide its rhs by d^/ds = 2n ■ p to get 

udF d£ 

d , (s) _ n djds C _ d f n»F(t(8))C \ 

ds P[) ~ 2n-p ~ds\ 2n-p J ' { <) 

where we omitted the argument s in 2n ■ p or d£/ ds since these quantities are constants 
of motion and used the fact that [n ■ p, £] = in 2 = 0. We also omitted the argument of 
dF/d^ and used that C is an s-independent operator. A direct integration of the above 
equation leads to 

^) = ^« + ^, (28) 
2n ■ p 

where is a constant operator that satisfies 

nrD , =ni ^ = m^m. (29) 

For simplicity, from now on, we shall omit indices and a matrix notation will be assumed, 
so that the previous equation is written simply as 

Two relations involving the constant operators C and D can be written from our previ- 
ous results. Inserting (28) into (26), and computing p^p^ from last equation we obtain, 
respectively, 

C = DaD (31) 
p 2 = D 2 + CF. (32) 

Hence, from equations (26) and (32), the proper-time hamiltonian (32) can be written as 

H = —D 2 . (33) 

Then, inserting (30) into (20) and integrating in s following a procedure analogous to that 
used to obtain p(s), we get 

, (s )-, ( 0H- 4 «M>-- 4 <«°>> fi^_ 2a0 W, (34) 
2n ■ p \n ■ p J 

where we defined A by 

Solving the previous equation for D, we obtain 

D _ M x(s) - x(0) _ nC A(Z(s))-Am) + 0(fl2) m 
D ~ M 2s 2n-p C(s)-m + i h { ) 

5 



where 



M :-- 



I — a 



Equations (31) and (36) allow us to write constant C in the form 



C= ~^(x(s) -i(0))aM 2 (i(s) -x(0)) + 0{a 2 ) . 



(37) 



(38) 



Using equations (33) and (36), and keeping terms only ut to order a, the proper-time 
hamiltonian in the weak field approximation takes the form 



H = -^(x(s)-x(0))M(x(s)-x(0)) . 



(39) 



To write the above proper-time hamiltonian in the appropriate s-ordered form we note 
the following commutation relations: 



K(0), = - 2n • p(s), x»(s)} = -2in»s , 



[M^x x (s),x li (0)] = 2s 



2n . P m-m 

Integrating equation id s {x's\x"0) = (x's\H\x"0), we get 

3>(x', x 



8is 



(40) 
(41) 



( x 's\x"0) = {X J 2 X ' exp [-^{x' - x")M{t', i"){x' - s")} , 



where 



I -a 



-1 



(42) 



(43) 



and x") is an s-independent quantity to be determined by imposing constraints (19). 
Let us impose the first constraint written in (19). To evaluate its rhs we need to compute 
d'(x's\x"0), while to evaluate its lhs we need the s-ordered matrix element of operator 
p(s). Differentiating equation (42), we obtain 



where 



.d'(x's\x"0) _ .<9'$ 
(x's\x"0) ~ 



d'M(?,e 



+ M 



x — x 



^ +-tf-3r)&M(?,g')tf-ar), (44) 



n 



no 



A(?) - A(C) 



aM 2 . 



(45) 



Using the solution for p(s), obtained from equations (30) and (31), a straightforward 
calculation yields 

d'®(x',x") = 0. (46) 



A similar calculation gives 



&'<$>{x',x") = 0. 



(47) 
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Hence $(x',x") is a constant, denoted simply by $. This constant may be evaluated by 
imposing the initial condition (18). However, the explicit calculation is very similar to 
that verifying the consistency condition in section 5, so that we avoid it here. The result 
is given by 

*=(SF- (48) 

Collecting all previous results, we finally obtain 



G{x', x" 



1 



ds 



lQn 2 Jo s 2 



exp ^ — -^-\ x ~ x ) 



I — a 



[X — X ) 



(49) 

As expected, if we take in the above expression the limit a — > 0, we recover the well known 
free particle propagator in the Minkowisk space. 



3 Proper-time method for a Dirac particle 

For a spin 1/2 particle the Green's function satisfies the equation 

- % -F{i')a^ ltl dl - m) S F (x',x") = 8(x' - x") , (50) 

where 7 M are the usual Dirac matrices. As it is common in problems involving fermion 
Green functions in external fields, we conveniently define Af(x',x") as follows: 

S F (x',x") := {i-fdl-^F^ar^dl + r^j A F (x',x"). (51) 
As a consequence, A F (x, y) satisfies a boson-like second order differential equation, namely, 

- FMardfll + -— <j^* v % + m 2 \ A F (x', x") = -6{x' - x") , (52) 

where = j v ] and we neglected terms quadratic in a. Using the notation 

introduced previously, A F (x',x") can be written as 

A F (x',x") = -i \ dse- m2s (x'\e- iHs \x") = -i / ds e~ im2s {x' s\x"$) , (53) 
Jo Jo 

where now the corresponding proper-time hamiltonian can be written in the form 

H = -p 2 + papF(0 + Ipacrn^- , (54) 

where we are following the same notation as before, namely, primed quantities are eigen- 
values while unprimed ones are operators. In writing last equation we also used the fact 
that \pap,F{£)\ = and [paan, dF(£) /d£] = 0. 

From now on, whenever it does not cause any confusion, we shall omit the argument 
s from the operators involved as well as the argument of F and its derivatives dF/ d£, and 
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d 2 F/d^ 2 . With this in mind, the Heisenberg equations of motion for the operators x and 
p are given, respectively, by 



dx 1 dF 

- = 2(1 - aF)p - -aan— , 

dp ( dF 1 d 2 F" 

is = n [ pap ^ + 2 paan ie 



In this case, the constants of motion are 



n ■ p = 



2s 

Ci = pap , 
C 2 = paan . 

Integrating the equation of motion for p we have 

/ \ n ( ^ 1 ^ dF s 

^ = *r- p { ClF+ 2 c ^ 



D 



(55) 
(56) 



(57) 

(58) 
(59) 



(60) 



where Df is a constant. For future convenience, we use the previous equation to write p 2 
in the form 

12 1 '"' ' "<"* " (61) 



j? = D) + dF + ^C 2 - 



1 rfF 

2^' 

where we used the fact that n-p — n- Df. Substituting equation (60) into the Heisenberg 
equation (55) and integrating in s, we obtain, after some convenient rearrangement, 



D f = 



A(gs))-A( m 



x(s)-x(0) \ l n _ F(tj(s))-F(m) 
2s J 4 £(*)-£(0) 



2n • p 



r ^))-^(Q)) i r *TOQ)-F(g(o)) 

, 1 i(s)-m + 2 2 e(«)-e(o) 



+ 0(a 2 ) . 



(62) 



Since n^a^ = 0, and keeping only terms up to first order in a, constants C\ and C 2 are 
given, respectively, by 



Ci = D f aD f 
C 2 = Dfaan 



/ x(s)-x(0)^ fx(s)-x(0)' 
x(s) — x(0) 



2s 



\ 2s 
aan + 0(a 2 ) . 



+ 0(a 2 ) 



(63) 
(64) 



From equations (54) and (61), the proper-time hamiltonian for the case of a Dirac particle 
is given by 

H = -D 2 f , (65) 
which, in the linear approximation (weak gravitational field), can be written as 



H 



x(s)-x(oy 

2s 



I + a 



A(£(s)) - A(Z(0)) 



x(s)-x(0Y 
2s 



(66) 



One may repeat the steps followed in the last section to get 



(x's\x"0) = 



exp I - — (x' -x") 



I+a M?)-A(n 



(x 1 - x") , (67) 



where <3>(x', x") is to be determined. The constraints to be satisfied are the same as in the 
previous case, namely, (x's\p(s)\x"0) = id'(x's\x"0) and (x's\p(0)\x"0) = —id"(x's\x"0). 
The first constraint leads to 



<9'$ 



d>-f^(F(Z>)-F(e)) 
An ■ p 



where 



(C 2 ) := (x's\C 2 \x"0) = ^' X "\ an + 0(a 2 ) . 



Hence, we can write 



§(x',x") = x(x")exp 



2s 



' 8{C 2 )(F(?) - F(?)) 



2(e-eo 



(68) 



(69) 



(70) 



where the ^''-dependence of x( x ") is determined by the second constraint written previ- 
ously. The imposition of this second constraint leads to the differential equatin d"x = 0, 
so that x is a constant, denoted by C : 



$(a/,x") = C exp 



' s(c 2 )(F(e)-F(n) 
2(e - e) 



(71) 



As before, we evaluate the remaining constant Cq by taking the limit s — > and the final 
result for A F (x',x") is written as 



A F (x',x") 



1 f°° ds _ im 2 s 

^ / e im s exp 

167T 2 Jo s 2 



4s 



[x — X 



I + a 



x exp 



(x' - x")aan{F{i') - F{£"))' 



(x'-x")^ 
(72) 



The desired fermionic Green function is then obtained by inserting last expression into 
equation (51). 



4 Green functions by operator techniques 

In this section we show that using operator techniques one can obtain the Green functions 
in a simple manner. We were motivated by successful calculations previously made for 
Green functions of relativistic charged particles in similar external electromagnetic fields 
[7, 8]. We shall start by discussing the case of a scalar particle. The essential idea is to start 
with expression (16) for the Green function, namely, G(x',x") = —i ds e~~ im s (x's\x"0) 
and then relate by an algebraic method the expression of Schrdinger-like propagator 
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(x's\x"0) for the problem containing the interaction with the corresponding one for the 
free case, denoted by (x's\x"0) . For the scalar particle we note that 



p 2 -papF(O = S p 2 S - 1 , 
where we defined the operator S by 



So = exp 



.papA(£) 



2n ■ p 

and used the well known relation for operators A and B 

e A Be- A = ^ C n , where C = B and C n+1 = [A, C n \ 

n 

As a consequence, we have 



(73) 



(74) 



(75) 



/\^is(p 2 -pa P F)\ x ir 



(x e 



(x'\S e lsp Sq-^x") 



(x e 



isp 



(e- isp2 S e l 



lsp Sq-'Ix 



= (x's\S (s)S (0)- l \x"0) , 



(76) 



where we conveniently inserted the identity operator e lsp e %sp and the meaning of the 
subscript o is such that 



(x's\x"0) o 



(Aits) 



Inserting (74) into (76), a straightforward calculation leads to the matrix element 



(x's\x"0) 



(Atts)' 



exp < (x — x 

' ! 4s v 



I + a 



A(£) - A(?') 



[x — X 



(77) 



(78) 



as an exact result. It is also the same as that found by the proper-time method in the 
small a approximation. This is curious. Perhaps the reason is that we have imposed 
very few restrictions on a. In the case of the electromagnetic field the tensor and its 
dual have special properties which allow simplifications. The other difference is that the 
proper-time equations of motion are nonlinear. 

For the Dirac particle we may proceed in an analogous way. With this purpose, we 
now define operator S as 

r i paan 



exp 



4 n ■ p 

Using the BCH-like relation (75) it is straightforward to show that 

i dF 

p 2 - papF(^) - ijPaon— = S(p 2 - papF^S' 1 . 



(79) 



(80) 



Note that operator S eliminates the interaction term containing explicitly the Dirac ma- 
trices. To reduce further the remaining term can be done with the aid of operator S , 
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introduced before for the spinless particle. Hence, we obtain for the Dirac case, the 
analogous expression of equation (76), namely, 



(x'\exp < is 



2 i dF 



\x") = (x's\S(s)So(s)(S(0)So(0))- 1 \x%)o 



(81) 

Inserting in the previous equation the expressions of S(s), S (s), S o (0) 1 and S(0) 1 it is 
not difficult to obtain the desired Green function for the Dirac particle. 

We close this section by showing an alternative way of employing albegraic methods 
in order to obtain the above relativistic Green functions. As we shall see, the present 
method leads directly to the Green functions without using the proper-time method. 

Starting with the spin case, we have 

G(x',x") = (x'\G\x") , (82) 

where 

G = (p 2 - pa P F(0 - m 2 y l = S (p 2 - m 2 )' 1 ^" 1 , (83) 
with operator S defined by equation (74). Hence, we can write 

tAA Y-So-'Hd") (84) 

27T 4 v — m l + ie 



which leads to the final result 



For the spin | the Green function Sf(x',x") is given by 

S F (x', x") = (i-fdl - l -F(?)a^dl + m) A F (x', x") , (86) 
where Af(x',x") is given by 

^f(x',x")= f f P f^'^. exp\-i^[A(a-A(n})expU^[F(x')-F(x'')}), 

(87) 

which agrees with the result obtained in section 2.2. 



5 Ansatz solution for the Green function 

In this section we presen an alternative procedure to construct previous Green functions 
which is based on a convenient ansatz for the desired solution. We ilustrate the method 
for the spin particle but a generalization to the spin 1/2 case can be made without 
difficulty. The main motivation for such a procedure is that it worked very well for the 
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caso of a relativistic charged particle under an external electromagnetic field of a plane 
wave [15]. The desired Green function satisfies the differential equation 

- a^F{i)d"d v + m 2 ] G(x - x') = -8(x - x') , (88) 

where, for convenience, we are using in this section variables x and x', instead of variables 
x' and x" chosen in the preceding section. Hence, along this section, x is not an operator. 
As usual, the Green function can be written as 

G(x, x') = r ds e~ rm2s A{x, x', s) , (89) 
Jo 

where A(x,x',s) satisfies the Schr dinger- like differential equation 

[id s - + a^F^Sr] A(x, x', s) = 0, (90) 
subjected to the initial condition 

iA(x,x',s) s ^ 5(x-x'). (91) 
Now, in order to factorize the free particle solution, we try the following ansatz 

A(x, x', s) = A (x, x', s)E(x, x', s) , (92) 
where A (x,x') corresponds to the free particle solution, that is, 

(id s -d 2 ) A (x,x') = 0, (93) 
whose the well known solution is given by 

A ° fel ' s) =16^ eXP {"^ ! }- <94> 
Substituing our ansatz (92) into equation (90), we have 

A i<9 s E - 2^A ^S - A <9 2 E + Fa^ (d^d u A ) E + 

+ 2Fa^A d u Z + Fa^ (^9"E) A = 0. (95) 

Since h^ u is small, we neglect the last two terms in above and get 

A i<9 s E - 2<9 M A <9 M E - A <9 2 E + Fa^ (d^d u A ) E = 0. (96) 

Next, using that 

d,A = -*^^Ao 
d,O u A = ^A - {X ~ X,) ^- X% A (97) 
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equation (96) takes the form 



E = 0. 



(98) 



An inspection of the previous equation suggests us to try a solution of the form 

'-if(S,S')(x-x') ll ar(x-x') v \ 



S = exp 



4s 



(99) 



In order to substitute last expression into equation (98), we need 

<a , s = _/'-^('-'V 



4s 2 



and 



df, 



n^(x - xXcTlx - x') u + 2fa^(x - x'f 



E. 



(100) 
(101) 



Further, we also have 



^" E = "16? 



n„— (x - x') a a a(S (x - x')p + 2fa^ a (x - x' 



x 



x 



n^(x - x') a a^{x - x% + 2faT(x - x') a 



where we have used the properties of and a^ u , namely, n 2 = and n^a^ u = 0. 
order of magnitude of our interest we may write 



d^H = 



Hence, substituing equations (103), (101) and (100) into equation (98), we get 

d 



We have then 



where we defined function A by 



dA 



(102) 
To the 

(103) 

(104) 

(105) 
(106) 



Therefore, equation (92) takes the form 



A(x, x', s) 



16"7T 2 S 

i 

167T 2 S' 



- exp < -i 



. (x - x'f + (x- xYa^Q^ix - x'f 



4s 



[X 



exp < — %- 



- x'f + (x - x'Ya ^llfP (x - x'Y 



4s 



(107) 
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With the purpose of checking the self-consistency of the result just obtained, not that 

-1^" A m \ ;(„\2 



lim / d A k e ^(B^s-A^) = Um \_. A Xll (B Ax u 1 t {nf 
■^ + J o + I 4s J s ^det(B) 

As a consequence, we have 



i^v exp {-" &x " {B ^r Ax " } - ^ ~ < 109 > 



Making the following identifications, 

AA 



B"„ = 5", - a"„^ 

we get 

det(B) = 1 + Tr(a\) A ^~p^ = 1 , (110) 
which confirms the initial condition in the parameter s, 

iA(x, x', s) s -=^ S(x - x') . (Ill) 

Inserting expression (107) into equation (89), we obtain the desired Green function, in 
agreement with our previous calculations. 

6 Semiclassical approximation 

As we have seen in many of the previous calculations, relativistic Green functions may 
be written in terms of Schrodinger-like propagators, if we introduce appropriately an in- 
tegration over the so called proper-time s. For instance, the Green function for a scalar 
particle can be written as Gf(x',x") = —i Jq°° dse~ ls( - m ~ le "> (x's\x"0) (see equation (7)), 
where (x's\x"0) can be interpreted as a Feynman propagator of an auxiliary problem of 
a non-relativistic particle in four dimensions whose dynamics corresponds to the evolu- 
tion in the parameter s, which plays the role of time in this auxiliary problem. Once 
(x's\x"0) behaves like a non-relativistic Feynman propagator, we have at our dispposal 
all techniques developed to compute this quantity as, for example, the Feynman path 
integral method. In particular, in the context of path integrals, it may be convenient to 
use the semiclassical approximation. It is well known that whenever the corresponding 
lagrangian is quadratic in the coordinates and velocities the semiclassical result gives the 
exact result. 

The purpose of this section is to show that if we apply the semiclassical approximation 
to the problem at hand (a relativistic particle in the background of a plane wave gravi- 
tational field) we shall obtain the exact result. At first sight this is an unexpected result, 
since the lagrangian of the corresponding classical problem is far from being quadratic. 



14 



However, there is a strong reason which suggests that this will be the case, namely, the 
semiclassical method when applied to the problem of a relativistic charged particle in an 
external field of an electromagnetic plane wave yields the exact result. Here, we shall dis- 
cuss only the case of a scalar particle. To avoid any confusion with the notation, observe 
that, in this section, all quantities are not operators, but classical numbers. 

The semiclassical approximation for the Feynman propagator (x's\x"0) is given by 



(x's\x"0) 



{2m)' 



d 2 S r 



dx'dx" 



1/2 



o iS cl 



where S c i means the functional action 

S(x,j) = L(x A1 (r),x M (r)) dr 



;i!2) 



(113) 



evaluated with the classical solution x^, that is, S c i = S(x c i), where x c i satisfies the 
Euler-Lagrange equations 



d ( dL 
ds \ dx„ 



dL 



dx, 



/i = 0, 1,2,3, 



submitted to the Feynman conditions 

*5(r = 0) = s w " ; xS(r = a) 



x 



in 



;ii4) 



(115) 



Hence, in order to use this approximation, we need to construct the classical lagrangian 
corresponding to the following classical hamiltonian 



H(x,p) = -p 2 +p /1 h^ u (x)p" 



(116) 



Recalling that p^ = dL/dx^, it is not difficult to show that the corresponding lagrangian 
can be written as 



1 1 

L(x,x) = --{x^ + x^h^x)^ = --g iiv [x)x i "x h 



The functional action is then given by 

1 

'o 



S(x) = /J [-li^(r)^(r) - - A x»{T) ailv x v {T)F{n ■ x{r)) 



dr 



(117) 



(118) 



The application of Euler-Lagrange equations (114) to the lagrangian (117) lead to 
classical equations completely analogous to the Heisenberg equations discussed in section 
2, but do not forget that here x^ and p^ ara not operators. However, exactly the same 
kind of solutions are obtained from these equations, so that we just write the classical 
solution as 



*50-) = s5(o) + 



2(n • pY 



n ■ p 



;il9) 
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where C and are (classical) constants to be determined by imposing Feynman condi- 
tions (115). From the classical version of (38) we have, up to first order in a, 



_ 1 (x' - - x") v _ 1 Ax^Ax" 

6 "4 ? " : 4 ? ' (120) 



where we defined Ax^ 1 = (x' — x")^. 

In order to obtain D^, we take r = s in equation (119) and contract x^(s) with n M 
and a^ v to obtain, respectively, 

n^Ax, = -2D»n u s =>- n - D — -— ; (121) 

a^Ax, = -1D v aTs, a^D v = — — , (122) 

where we defined A£ = n^{x' — x")^ — £' — Substituting these relations into the 
solution (119) with r = s, we obtain in terms of x' and x", namely, 

where we AA := A(£') — A(£"). Substituting expressions (120) and (123) into equation 
(119), and then differentiating with respect to r, we obtain the classical velocity x^(r) at 
any "instant" r in terms of x' and x", namely, 

Substituting last equation into the functional action (118) and keeping only terms ur 
to first order in a, we obtain after a lengthy but straightforward calculation the desired 
classical action, 

Now we need to evaluate the Van Vleck-Pauli-Morette determinant to obtain the pre- 
exponencial factor of the semiclassical propagator. With this purpose, first note that 



d 2 S r 



dx' fix" v 2s 

As a consequence, we have 



UfT + Oia)). (125) 



det (a^) = is* det (" + ° (tt » = ^ 1 + TYO(tt) ) ■ (126) 
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However, note that tr(0{a)) = 0{a 2 ) so that det 1/2 (^§^7) = 4^2 + 0{a 2 ). Collecting 
the previous results and using equation (112) we conclude that, up to first order in a, the 
Feynman propagator (x's\x"0) in the semiclassical aproximation is given by 

Substituting last expression into equation (16), we reobtain the correct Green function 
for a scalar relativistic particle in the weak gravitational field of a plane wave. 



7 Conclusions and final remarks 

In this work we have presented a few alternative techniques to calculate Green functions 
of relativistic particles under the influence of a weak gravitational field of a plane wave, 
with particular attention to the Fock-Schwinger proper-time method. In fact, we started 
this paper by a detailed application of this method in the construction of the Green 
functions for both spin zero as well as spin 1/2 particles. We showed how these Green 
functions can be obtained by algebraic methods in an extremely compact and elegant 
way. We also showed how to reobtain these solutions with an appropriate ansatz for 
the Green function and, finally, we discussed a semiclassical solution and checked that 
for the case at hand this approximation is sufficient to give the exact result. From one 
hand, this is a surprising result, since the lagrangian involved in the solution is far from 
being quadratic. On the other hand, it could have been guessed, for a semiclassical 
approximation yields the exact Green function for a relativistic charged particle in an 
external electromagnetic field of a plane wave. The results obtained here corroborates and 
complements those of reference [13]. We think that our calculations may be useful in other 
contexts. Particularly, the search for other non-quadratic problems whose exact solutions 
coincide with those obtained by the semiclassical approximation is and interesting issue 
and can tell us more about the role of such an approximation in physical problems. 
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